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1. Introduction 

For the past ten years or so, a rumour has been making the rounds of particle physics 
Conferences and Workshops: measured bottom quark production rates are significantly 
larger than predicted by next-to-leading order QCD. 

Such a state of affairs casts doubts on the ability of QCD to properly describe this 
process, and opens the way to speculations that either drastic improvements to QCD 
calculations are needed, or even that effects beyond the Standard Model might be at work. 

We shall challenge the rumour under two main aspects, taking aim especially at the 
words "measured" and "significantly". It is our hope that, once done, we shall have con- 
vinced the reader that standard approaches to QCD calculations do a good job in describing 
the experimental findings. 

Let us first briefly review the history of the experimental results. Measurements of 
the bottom transverse momentum [px) spectrum at colliders began in the late 80's, when 
the UAl Collaboration, taking data at the CERN SppS with = 546 and 630 GeV, 
published results for the px > rrib (the bottom quark mass) region [1,2]. These results 
were compared to the then recently completed next-to-leading order (NLO), i.e. order af, 
calculation [3-6], and were found to be in good agreement. 

During the 90's the CDF [7-13] and DO [14-17] Collaborations also measured the 
bottom quark pT distribution in pp collisions at the Fermilab Tevatron at y/S = 1800 GeV. 
Apparently at odds with the UAl results, the Tevatron data seemed to display an excess 
with respect to NLO QCD predictions. At the same time, rates for bottom production that 
appeared higher than QCD predictions were also observed in 77 collisions by three [18-20] 
LEP experiments, and by the HI [21] and ZEUS [22] Collaborations in ep collisions at 



HERA. 
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By the end of the millennium the 'excess' was therefore apparently so firmly established 
that the experimental papers showed no shyness in proclaiming it to the world, e.g. "T/ie 
differential cross section is measured to be 2.9±0.2 (stat©systpj,)ib0.4 (systj^) times higher 
than the NLO QCD predictions..." [12] and "...NLO QCD calculations underestimate b 
quark production by a factor of four in the forward rapidity region" [16]. 

Despite this seemingly overwhelming evidence, we shall argue that QCD is instead 
rather successful in predicting bottom production rates. Improved theoretical analy- 
ses [23, 24] and more recent experimental measurements by the CDF [25] and ZEUS [26] 
Collaborations support this claim, which is also borne out by a critical reconsideration of 
previous results. 

2. The Paradigm 

We shall take NLO QCD calculations as a benchmark for comparisons. We shall require 
the experimental measurements to be genuine observable quantities. By this we mean that 
as a matter of principle we do not wish to compare 'data' for, e.g., 6-quark px distributions, 
since such a quantity is clearly an unphysical one: the quark not being directly observed, 
its cross sections have to be inferred rather than directly measured. 

A meaningful comparison will therefore be one between a physical cross section and 
a QCD calculation with at least NLO accuracy. Non-perturbative information, where 
needed, will have to be introduced in a minimal and self-consistent way. This means 
that we refrain from using unjustified models, and we shall only include non-perturbative 
information that has been extracted from one experiment and then employed in predicting 
another observable, using the same underlying perturbative framework in both cases. Such 
a precaution allows for a good matching between the perturbative and the non-perturbative 
phases, a necessity in that only the combination of the two steps leads to an unambiguous 
measurable quantity. 

In practice, the non-perturbative information relative to the hadronization of the b- 
quarks into S-hadrons is extracted from LEP data with a calculation [27] which has NLO -|- 

NLL accuracy^. The framework presented in [23] is used: the LEP (or SLD) data [28-31] are 
translated to Mellin moments space, and only the moments around N = 5 are fitted. This 
ensures that it is the relevant part of the non-perturbative information which is properly 
determined.'^ These non-perturbative moments are then used together with a calculation 

'^By NLL accuracy we mean that large terms of quasi-collinear origin, proportional to log" (Q^/m^) 
and a" \og"~^{Q^ /ml), Q being the LEP centre-of-mass energy ~ 91.2 GeV, are resummed to all orders. 

^One readily realizes [32,33] that these moments arc the only important ones when convoluting a trans- 
verse momentum spectrum which falls off like dai/dpT — A/p^, as is the case in hadronic collisions. We 
have in fact that the convolution of such a perturbative spectrum with a non-perturbative fragmentation 
function gives 

das d&b ^ /X , [ dz { z \ '' ^ , , dijb „r 
apT dpT J z \pt J dpT 

i.e. the spectrum for B hadrons is given by the spectrum for 6 quarks multiplied by the fifth moment of 
the non-perturbative fragmentation function. 
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having the same perturbativc features, FONLL [34] (Fixed Order plus Next-to-Leading 
Log - in this case log(p^/m^)), to evaluate the cross sections in pp collisions. 

The expectation is then that total cross sections be reproduced by the NLO calculations 
for b quarks, and that differential distributions for B hadrons be correctly described by 
a proper convolution of the FONLL perturbative spectrum for b quarks and the non- 
perturbative information extracted from LEP data. Notice that a minimalist use of non- 
perturbative information is made: there is no attempt to fully describe the hadronization 
process. Only the relevant phenomenological information is determined from data and used 
in the predictions. 

A successful comparison will see data and theory in agreement within their combined 
uncertainties. The theoretical ones will be assessed by varying as extensively as reasonable 
the parameters and the unphysical scales entering the predictions. As for the experimental 
errors, it is perhaps worth reminding that only 1-sigma errors are usually shown on the 
plots, so that non-overlapping bands do not necessarily point to a solid disagreement. 

3. The Data 

3.1 77 

Let us first consider the gamma-gamma data [18-20], measured in e+e^ collisions at = 
194 GeV. Unfortunately only the L3 paper is published in final form. The results from 
the three experiments appear fully compatible, reading a{e^e~ e'^e~bbX) = 12.8 it 
1.7 ± 2.3 pb (L3), 14.2 ± 2.5 'fl pb (OPAL) and 14.9 ± 3.3 ± 3.4 pb (DELPHI). These 
results should be compared to a theory prediction that the experimental papers estimate 
of the order of 4 pb, apparently with an uncertainty inferior to 10-15%. Figure 1 shows 
the comparison between theory and experiment as usually presented at Conferences. 

Four comments are in order here. The first is that the theoretical uncertainty has 
probably not been fully explored, missing at least the investigation of the effect of different 
photon parton distribution function sets, which are only approximately known, and the 
exploration of independent renormalization and factorization scale variations. Ref. [35] 
for instance, while not performing a full exploration of all variables, still estimates a larger 
uncertainty. The second is that, especially when including theoretical errors, the real 
distance between data and theory probably does not exceed 2 sigmas. Hence, claims of 
data being "...in excess of the QCD prediction by a factor of three" [18] appear premature, 
not being backed by adequate significance. An 'excess factor' may correctly refer to central 
values, but it should not be used to measure a 'discrepancy' until errors are also included. 
The third comment is that the three experiments essentially relied on the same technique 
and the same tool (the PYTHIA Monte Carlo) for extracting the b signal. Hence, their 
results are probably strongly correlated and their accord is probably less significant than it 
appears. Finally, last but not least, the total cross section is not what is actually measured 
by the experiments. Rather, only a fairly tiny fraction of the cross section is measured 
with the use of PYTHIA, and by means of the same Monte Carlo this cross section is then 
extrapolated to the full phase space. Needless to say, such a procedure might introduce a 
bias, as the theoretical prediction of PYTHIA outside the measured region cannot of course 
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Vs (GeV) 

Figure 1: Total cross sections for 77 cc and bb. 

be considered a priori correct. Unless the uncertainty that the extrapolation introduces 
can be reliably estimated, the total cross section results cannot be considered as real 
'measurements'. Comparisons of these data with theoretical predictions should therefore 
be taken with a grain of salt. 

3.2 7p and ep 

Bottom production cross sections have also been measured in both photoproduction and 
Deep Inelastic Scattering (DIS) at HERA by the HI and ZEUS Collaborations. The first 
measurement was performed by the HI Collaboration [21] in 1999. While allowing for 
experimental and theoretical uncertainties of the order of 10-20%, this paper (and in par- 
ticular the erratum) presented a cross section (deconvoluted and extrapolated) about a 
factor of three larger than the NLO prediction. Shortly thereafter the ZEUS Collaboration 
also published its first data on bottom photoproduction [22]. Comparison with NLO QCD 
was again performed after deconvolution to parton level and extrapolation by Monte Carlo 
to a more inclusive cross section. While the central value of the cross section was found to 
be about a factor of two larger than the central NLO prediction, large experimental and 
theoretical uncertainties should prevent one from inferring an incompatibility between the 
two numbers, which should actually be considered compatible at the 1-sigma level. Never- 
theless, the paper concludes with the observation that the result is ...consistent with the 
general observation that NLO QCD calculations underestimates beauty production..." . 

The reliability of these two results being somewhat diminished by the deconvolutions 
and/or extrapolations performed by Monte Carlo before actual comparisons to NLO QCD 
predictions, the ZEUS Collaboration set out to perform direct comparisons between theo- 
retical predictions and real measurements. Photoproduction data [26] are found to be in 
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Figure 2: One of the plots published in [26] of the observable cross sections measured by ZEUS 
and their comparisons to theoretical predictions. 



very good agreement with NLO QCD (see figure 2), prompting the conclusion that the large 
excess found by HI is not confirmed. It is worth noting that the ZEUS Collaboration also 
mentions that the new results are consistent with their old one [22] . While this is certainly 
true from the statistical point of view, it becomes then apparent that different conclusions 
('above theory' vs. 'consistent with theory') were drawn from compatible measurements. 
One cannot help noticing that - apparently - a different weight is given to 'sigmas' when 
comparing experiment to experiment rather than experiment to theory. 

The ZEUS Collaboration has very recently extended the photoproduction analysis to 
the DIS regime [36]. The data are found to be consistent with NLO QCD predictions, and 
at most two sigmas higher in a few bins of the distributions. 

3.3 pp 

Hadronic collisions were historically the first to produce bottom production data in collider 
physics. The first results date back to 1986, and were published by the UAl Collaboration, 
using pp collisions at ^/S = 546 and 630 GeV. In a subsequent paper [2], the UAl Collab- 
oration compared bottom quark transverse momentum distributions to the then recently 
completed NLO prediction [3-6], and found a good agreement (see figure 3a). It is worth 
noting that, while comparing data deconvoluted to the unphysical quark level, the paper 
also included data for the observed B mesons. While a comparison with theoretical predic- 
tions at the hadron level was not performed at that time, it is however still possible now, 
thanks to the preservation and publication of the real data. Such a comparison, performed 
with the modern tools and framework described in Section 2, shows an agreement similar 
to the one observed at the quark level. 

The main source of bottom data in pp collisions in recent years has been the Fermilab 
Tevatron, running at VS = 1800 GeV first and 1960 GeV later. The first results from 
this machine were given in 1992 by the CDF Collaboration [7]. The inclusive cross section 
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(integrated above a minimum of 11.5 GeV, and within the central rapidity region \y\ < 1) 
was published deconvoluted to the quark level. Its central value was found to be a factor 
of six larger than the NLO prediction, but the very large errors only made it a 1.6 sigma 
distance and therefore not a significant one. 

One year later CDF started publishing [8, 9] the plot of the transverse momentum 
distribution which will then become the icon of the supposed 'excess' (see figure 3b). The 
data were published only at the unphysical quark level and, while reporting differences 
of the order of 1-2 sigmas between data and theory, the papers still conclude that "...the 
next-to-leading order QCD calculation tends to underestimate the inclusive 6-quark cross 
section" , hinting therefore for the first time that a discrepancy might be present. The same 
conclusion was reached in a subsequent paper [11], where cross sections for real particles, 
B mesons, were finally published. In this case the 'disagreement' was quantified for the 
first time, a fit to an overall scale factor between data and theory yielding 1.9 it 0.2 it 0.2. 

Around the same time CDF published the first data on bottom quarks, the DO Collab- 
oration also released some preliminary results which were presented at a number of confer- 
ences. Somewhat at odds with the CDF ones, they were in very good agreement with the 
NLO predictions. Therefore, weighing the data from both collaborations, speakers at the 
conferences (see e.g. [37,38]) generally reported a good agreement between bottom quark 
data and NLO QCD. 

Given this state of affairs, the final DO data must have caused some surprise when, 
published in final form about one year later [14], they became more CDF-like, now lying 
around the upper edge of the theoretical uncertainty band. The prediction of NLO QCD, 
however, was still considered to be giving an adequate description of the data. 

The subsequent pair of DO papers on this subject [16, 17] should have caused an even 
larger surprise. Despite the conclusions of the previous paper ("adequate description" [14]), 
in the Introduction of [16] the previously measured b quark cross section is now considered 
to have been found "systematically larger" than the central value of NLO QCD predictions. 
This helps to digest the news that the data now show a considerable excess: "T/ie ratio 
of the data to the central NLO QCD prediction is approxim,ately three..." . [16] This strong 
statement is even upped in the following paper which, already in the Abstract, states that 
We find that next-to-leading- order QCD calculations underestimate h- quark production by 
a factor of 4 in the forward rapidity region" . [17] 

If we are to take such statements^ at their face value, we are of course to conclude 
that by the year 2000 the 'excess' appeared firmly established. Even new physics, in 
the form of light supersymmetric particles [39], was advocated in order to try explaining 
the discrepancy (see however also [40], which excludes this model using e^e~ hadronic 
data). More conventional explanations involve trying to consider small-x resummation 
effects [41]. These results hint that important contributions might be present, and certainly 
deserve further investigations. However, for the time being they lack overall NLO accuracy, 

^In passing, we note that no uncertainties are explicitly included, making these statements somewhat 
void of significance. Where errors are considered, often buried deep in the text and tables in the papers, 
they usually neglect the theoretical uncertainties, leading therefore to an overestimation of the significance 
of a possible gap. 
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Figure 3: A collection of bottom quark production measurements and comparisons with theory. 
From top to bottom, left to right: a) the UAl results [2], b) a collection of CDF and DO data, c) 
the CDF Run I i?+ cross section [12], and d) the DO &-jets cross section [42]. 



suffer from large normalization uncertainties, and do not satisfy the requirements of our 
'paradigm'. We shall therefore refrain from performing detailed comparisons using these 
results. 

Two further experimental papers were then added to the field: both are worth examin- 
ing closely, albeit for different reasons. The DO Collaboration performed a measurement [42] 
of the transverse energy distribution of jets containing a bottom quarks, the so called b- 
jets. These objects are real observables, while at the same time being largely independent 
of the fragmentation properties of the bottom quarks. By contrast, other observables like 
the B mesons do instead of course depend on the way the bottom quarks fragment into 
the bottomed hadrons. DO found that the cross section for 6-jets was compatible with the 
NLO QCD prediction [43] (see figure 3d). The CDF Collaboration updated instead [12] 
its measurement for the B mesons transverse momentum distribution, superseding [11] 
and representing the final analysis for bottom production with the Tevatron Run I data: 
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Figure 4: On the left (a), 5+ mesons production at the Tevatron and the new theoretical pre- 
diction. On the right (b), description of moments of B fragmentation data at LEP with different 
non-perturbative fragmentation functions. Both plots from [23]. 



"T/ie differential cross section is measured to be 2.9 ± 0.2 (stat © syst^^) ± 0.4 (systj^) 
times higher than NLO QCD predictions..." (figure 3c). A couple of comments are wortli 
making. The first is that, once more, the errors on this ratio do not include the theoretical 
uncertainty, which from perturbative sources alone would be at least of the order of 20-30% 
of the NLO result. The second is that the 'NLO QCD' prediction must of course include 
the non-perturbative information needed to fragment the bottom quark into the B meson. 
In this experimental paper this fragmentation was performed using the Peterson et al. [44] 
functional form, with its free parameter set to = 0.006, a standard choice dating back 
to determinations performed in 1987 by Chrin [45]. Such a procedure however neglects the 
notion that neither the bottom quark nor its fragmentation into B hadrons are physical 
observables. Neither of them is separately measurable, only their final combination is. It 
is therefore wrong in principle (and also, as we will see, in practice) to rely on a fixed and 
standard determination of the non-perturbative fragmentation function and to convolute 
it with whatever perturbative calculation is being used. The non-perturbative fragmenta- 
tion must rather be determined from data (usually from e~^e~ collisions) using the same 
perturbative framework and parameters (bottom mass, strong coupling) of the calculation 
which will then be employed to calculate bottom quark production in pp collisions. 

The problem of using the 'wrong' non-perturbative fragmentation function can become 
irremediable if only the data for the deconvoluted data for the unphysical b quarks are 
finally published. These data might in fact be biased by the deconvolution, but it would 
be very hard to reconstruct the originally measured ones. On the other hand, if the real 
measurements are (also) published {B mesons, 6-jets, muons or J/^''s from B mesons 
decays), they can be directly compared with theoretical predictions which include also the 
non-perturbative and decay stages. 

This is precisely what was done in [23]. The non-perturbative information was de- 
termined, as described above, from LEP data in moment space, and employed - within 
a consistent perturbative framework - to predict the B mesons transverse momentum 
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distribution at the Tevatron. Figure 4a shows how the data from [12] are compatible 
with the updated theoretical prediction: the ratio can now be roughly estimated to be 
1.7 ± 0.5 (expt) lb 0.5 (th). There is therefore no significant discrepancy between the data 
and the theory. Figure 4b shows how using the Peterson et al. form with = 0.006 does 
indeed underestimate the moments around N = 5, consequently leading to an underesti- 
mation of the rate in pp collisions. 

4. The Recent Comparisons 

The understanding of the potentially large biases related to the description of the non- 
perturbative fragmentation phase, and the proper inclusion of uncertainties from all the 
sources, allowed one to conclude that no significant discrepancies were probably present 
in the bottom data collected at the Tevatron during the Run I. These data were however 
always above a minimum pT of about 5 GeV. Since a harder or softer non-perturbative 
fragmentation function will leave the total cross section unchanged while shifting contri- 
butions to larger or smaller pT values respectively, it was still possible that such a shift 
was only faking a larger cross section. Small-p^^ data, and possibly a total cross section 
measurement, are therefore crucial for establishing whether the NLO QCD prediction does 
indeed account (or not) for the number of bottom quarks produced at the Tevatron. 

Such data, from the Tevatron Run II, have been recently made public in preliminary 
form by the CDF Collaboration [25] , and promptly compared [24] to the predictions given 
by the framework put forward in [23]. The data are in the form of J ftp's coming from bot- 
tomed hadrons Hi,. The theoretical predictions depend solely on the following calculations 
and parameters: 

• Perturbative inputs 

— FONLL calculation (i.e. full massive NLO calculation plus matching to NLL 
resummation) , both for e+e~ [27] and for pp [34] collisions 

— bottom quark pole mass = 4.75 GeV (varied between 4.5 and 5 GeV) 

— strong coupling (A^^) = 0.226 GeV, i.e. as{Mz) = 0.118) 

— renormalization and factorization scales (varied between fio/2 < ^lr^f < 2/xo, 



• Non-perturbative/phenomenological inputs 

— gluon and light quarks PDFs (CTEQ6M [46] default choice, MRST [47] and 
Alekhin [48] sets also used) 

— h quark to Hi^ hadron fragmentation (fitted to moments of LEP data, see [23,24]) 

— Hf, to J /ip branching ratio, 1.15% [49] and decay spectrum (from CLEO [50] or 
BaBar [51] Collaborations) 

After extensive exploration of all the numerically meaningful uncertainties, the predictions 
compare to the measured total cross sections as follows: 
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Figure 5: CDF J/tp spectrum from Hb decays, compared to theoretical predictions [24]. 
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The first two lines refer to physical cross sections, measured (and predicted) in the given 
visible region. The third line represents the deconvolution to the quark level. These results 
clearly indicate full consistency between theory and experiment within the uncertainties. 
The transverse momentum spectrum of the J/ip^s from 6's, shown in figure 5, is equally 
well described. The reason why the agreement now looks better than it did in figure 4 is 
twofold. On one hand, the theoretical prediction is increased by 10-20% by employing a 
more modern PDF set (CTEQ6M vs. CTEQ5M). On the other hand, the experimental 
data, which should have been about 10% higher due to the larger Run II energy, are instead 
about 25% lower (but still compatible with the old ones within the uncertainties). 



5. Conclusions 



Next-to-leading order QCD appears to be doing a good job in predicting total and single 
inclusive bottom quark production cross sections at the Tevatron and HERA. Comparisons 
performed at the observed hadron level, rather than at the unphysical quark level, do not 
seem to show significant discrepancies. Tevatron Run II preliminary results are even in 
very good agreement. We argue that discrepancies pointed out in the past were either not 
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very significant, in that the real size of the uncertainties might have been underestimated 
or simply overlooked, or that the 'data' might have been tainted by excessive use of Monte 
Carlo simulation in their extraction, deconvolution to parton level (perhaps with the wrong 
fragmentation function) , and extrapolation to full phase space. An extension of the number 
of comparisons performed at the observed particle level and in the visible regions will 
certainly help shed light on this point. 

New physics has been advocated at some point in order to explain the presumed 
discrepancy. While there is of course still room for it within the uncertainties, at the level 
of about 30% in the case of the Tevatron data, we argue that its presence is not needed in 
order to explain the single inclusive bottom production data. 

Finally, we wish to point out that much of the progress done in the last couple of 
years has been permitted by the possibility of comparing theoretical predictions to real 
data, rather than to deconvoluted/extrapolated ones. This is not always the case, as 
sometimes the original data are not published and are lost forever. We invite therefore 
the experimental Collaborations to always publish also results which do not depend (or 
depend as little as possible) on theoretical prejudices (e.g. in the form of a Monte Carlo 
code) for their extraction. This will avoid the risk of biasing them, and will leave open the 
possibility of performing updated comparisons in the future. 

Acknowledgment s 

Many thanks go to my collaborators, to Giinther Dissertori and Gavin Salam for the 
careful reading of the manuscript, to the experimentalists I talked to trying to understand 
the details of experimental measurements performed even years in the past, to all the 
Organizers for the invitations, and to those who raised their hands after my presentations at 
the Conferences, pointing out what papers I had overlooked or misinterpreted, and offering 
their opinion and advice. I have probably not convinced everybody with my arguments, 
nor perhaps retained all suggestions, but I have certainly enjoyed every conversation. 

References 

[1] C. Albajar et al. [UAl Collaboration], Phys. Lett. B 186 (1987) 237. 

[2] C. Albajar et al. [UAl Collaboration], Phys. Lett. B 256 (1991) 121 [Erratum-ibid. B 262 

(1991) 497]. 

[3] P. Nason, S. Dawson and R. K. Ellis, Nucl. Phys. B 303 (1988) 607. 

[4] P. Nason, S. Dawson and R. K. Ellis, Nucl. Phys. B 327 (1989) 49 [Erratum-ibid. B 335 

(1990) 260]. 

[5] W. Bccnakkcr, H. Kuijf, W. L. van Nccrven and J. Smith, Phys. Rev. D 40 (1989) 54. 

[6] W. Beenakker, W. L. van Neerven, R. Meng, G. A. Schuler and J. Smith, Nucl. Phys. B 351 

(1991) 507. 

[7] F. Abe et al. [CDF Collaboration], Phys. Rev. Lett. 68 (1992) 3403. 
[8] F. Abe et al. [CDF Collaboration], Phys. Rev. Lett. 71 (1993) 500. 



- 11 - 



[9] F. Abe et al. [CDF Collaboration], Phys. Rev. Lett. 71 (1993) 2396. 
[10] F. Abe et al. [CDF Collaboration], Phys. Rev. D 50 (1994) 4252. 

[11] F. Abe et al. [CDF Collaboration], Phys. Rev. Lett. 75 (1995) 1451 [arXiv:hep-ex/9503013]. 

[12] D. Acosta et al. [CDF Collaboration], Phys. Rev. D 65 (2002) 052005 

[arXiv:hcp-ph/0111359]. 

[13] D. Acosta et al. [CDF Collaboration], Phys. Rev. D 66 (2002) 032002 [arXiv:hep-ex/0206019]. 

[14] S. Abachi et al, [DO Collaboration], Phys. Rev. Lett. 74 (1995) 3548; 

[15] S. Abachi et al. [DO Collaboration], Phys. Lett. B 370 (1996) 239. 

[16] B. Abbott et al, [DO Collaboration], Phys. Rev. Lett. 84 (2000) 5478 [arXiv:hep-ex/9907029]. 

[17] B. Abbott et al [DO Collaboration], Phys. Lett. B 487 (2000) 264 [arXiv:hep-ex/9905024]. 

[18] M. Acciarri et al [L3 Collaboration], Phys. Lett. B 503 (2001) 10 [arXiv:hep-ex/0011070]. 

[19] OPAL Collaboration, OPAL-Note PN 455 (2000) 

[20] DELPHI Collaboration, Note 2003-015-CONF-635 

[http://delphiwww.cern.ch/pubxx/conferences/summer03/paper_aa0325.ps.gz] 

[21] C. AdlofF et al [HI Collaboration], Phys. Lett. B 467 (1999) 156 [Erratum-ibid. B 518 

(2001) 331] [arXiv:hcp-cx/9909029]. 

[22] J. Breitwcg et al [ZEUS CoUaboration] , Eur. Phys. J. C 18 (2001) 625 
[arXiv:hep-ex/0011081]. 

[23] M. Cacciari and P. Nason, Phys. Rev. Lett. 89 (2002) 122003 [arXiv:hep-ph/0204025]. 

[24] M. Cacciari, S. Frixione, M. L. Mangano, P. Nason and G. Ridolfi, JHEP 0407 (2004) 033 
[arXiv:hep-ph/0312132]. 

[25] CDF Note 6285; C. Chen, [CDF Collaboration], presentation at "Beauty 2003", 

http://www-cdf.fnal.gov/physics/new/bottom/030904.blessed-bxsec-jpsi/ ; M. Bishai [CDF 
Collaboration], presentation at Fermilab, Dec 5. 2003, 
http: / /www-cdf.fnal.gov/~bishai /papers / wandc.pdf 

[26] S. Chckanov et al [ZEUS Collaboration], arXiv:hep-ex/0312057. 

[27] B. Mclc and P. Nason, Nucl. Phys. B 361 (1991) 626. 

[28] A. Heister et al [ALEPH Collaboration], Phys. Lett. B 512 (2001) 30 
[arXiv:hep-ex/0106051]. 

[29] K. Abe et al [SLD Collaboration], Phys. Rev. D 65 (2002) 092006 [Erratum-ibid. D 66 

(2002) 079905] [arXiv:hep-ex/0202031]. 

[30] G. Abbiendi et al [OPAL Collaboration], arXiv:hep-ex/0210031. 

[31] G. Barker et al [DELPHI Collaboration], DELPHI 2002-069 CONF 603, 
http://delphiwww.cern.ch/pubxx/delnote/public/2002j069_conf_603.ps.gz 

[32] S. Frixione, M. L. Mangano, P. Nason and G. Ridolfi, Adv. Ser. Direct. High Energy Phys. 

15 (1998) 609 [arXiv:hep-ph/9702287]. 

[33] P. Nason et al, Report of the 1999 CERN Workshop on SM physics (and more) and the 
LHC, arXiv:hep-ph/0003142. 



- 12 - 



[34] M. Cacciari, M. Greco and P. Nason, JHEP 9805 (1998) 007 [arXiv:hep-ph/9803400]. 

[35] S. Frixione, M. Kramer and E. Laenen, J. Phys. G 26 (2000) 723 [arXiv:hep-ph/0002112]. 

[36] S. Chekanov et al. [ZEUS Collaboration], arXiv:hep-ex/0405069. 

[37] A. Zicminski [DO Collaboration], 

littp://www.slac. Stanford. edu/spires/find/liep/www?irn=4294769. Prepared for 8th Les 
Rencontres de Physique de la Vallee d'Aoste: Results and Perspectives in Particle Physics, La 
Thuile, Raly, 6-12 Mar 1994 

[38] K. A. Bazizi [CDF and DO Collaborations], FERMILAB-CONF-94-300-E, 

littp://www.slac. Stanford. cdu/spircs/find/hep/www?r=fcrmilab-conf-94-300-c, Presented at 
29th Rencontres de Moriond: QCD and High Energy Hadronic Interactions, Meribel les 
Allues, France, 19-26 Mar 1994 

[39] E. L. Berger, B. W. Harris, D. E. Kaplan, Z. Sullivan, T. M. P. Tait and C. E. M. Wagner, 
Phys. Rev. Lett. 86 (2001) 4231 [arXiv:liep-ph/0012001]. 

[40] P. Janot, arXiv:liep-ph/0403157. 

[41] H. Jung, Mod. Phys. Lett. A 19 (2004) 1 [arXiv:hep-ph/0311249]. 

[42] B. Abbott et al. [DO Collaboration], Phys. Rev. Lett. 85 (2000) 5068 [arXiv:hep-ex/0008021]. 
[43] S. Frixione and M. L. Mangano, Nucl. Phys. B 483 (1997) 321 [arXiv:hep-ph/9605270]. 
[44] C. Peterson, D. Schlatter, 1. Schmitt and P. M. Zerwas, Phys. Rev. D 27 (1983) 105. 
[45] J. Chrin, Z. Phys. C 36 (1987) 163. 

[46] J. Pumplin et al., JHEP 0207 (2002) 012 [arXiv:hep-ph/0201195]. 

[47] A. D. Martin, R. G. Roberts, W. J. Stirling and R. S. Thorne, Eur. Phys. J. C 28 (2003) 455 
[arXiv:hep-ph/0211080]. 

[48] S. Alekhin, Phys. Rev. D68 (2003) 014002 [arXiv:hep-ph/0211096]. 

[49] K. Hagiwara et al. [Particle Data Group Collaboration], Phys. Rev. D 66 (2002) 010001. 

[50] S. Anderson et al. [CLEG Collaboration], arXiv:hep-ex/0207059. 

[51] B. Aubert et al. [BaBar Collaboration], Phys. Rev. D 67 (2003) 032002 
[arXiv:hep-ex/0207097] . 



-13- 



